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FOREWORD 


Improvement  of  land  mobility  technology  — a major  goal  of 
the  U.S.  Army  Tank-Automotive  Command  (TACOM)  — requires  simula- 
tion of  the  interaction  of  the  terrain,  the  vehicle,  and  man. 

Under  TACOM  contracts,  Grumman  has  developed  mathematical  models 
of  rigid  wheel-soil  and  pneumatic  tire-soil  interaction  conceived 
as  alternate  submodels  in  the  AMC  (Army  Material  Command)  Mobility 
Model.  In  these  first  generation  mathematical  models,  the  terrain 
was  assumed  homogeneous  throughout  its  depth  and  characterized  by 
one  set  of  Coulomb  strength  parameters.  A great  number  of  tire 
tests  have  been  performed  in  homogeneous  soil  beds  at  the  U.S.  Army 
Engineers  Waterways  Experiment  Station  and  elsewhere.  These  have 
made  validating  these  models  possible  for  a great  variety  of  homo- 
geneous soil  conditions. 

In  the  field,  however,  soil  conditions  are  rarely  homogeneous. 
Precipitation  and  evaporation  give  rise  to  a moisture  gradient 
throughout  the  depth  profile  of  the  terrain  and  to  corresponding 
changes  in  the  strength  properties  of  soil.  Characterization  of 
such  soil  conditions  requires  the  consideration  of  Coulomb  strength 
parameters  that  vary  with  depth,  either  discretely  or  continuously. 
Theoretical  research  performed  under  the  present  contract  has  shown 
that  the  plasticity  theory  can  be  applied  for  the  determination  of 
stress  states  in  nonhomogene ous  soil  conditions.  In  this  report, 
the  results  of  this  research  are  presented  together  with  their  ap- 
plication in  the  tire-soil  model  for  the  consideration  of  nonuni- 
form and  layered  soil  conditions. 


iii 


ABSTRACT 


Mathematical  models  of  tire-soil  interaction  have  been  de- 
veloped for  nonhomogeneous  soil  conditions  where  the  soil  strength 
varies  either  continuously  or  discretely  with  depth.  New  methods 
of  solving  the  differential  equations  of  plasticity  for  soils  have 
been  developed  for  the  bearing  capacity  problem  in  two-layer  soils. 
Composite  slip  line  fields  obtained  by  these  methods  and  the  asso- 
ciated bearing  stresses  are  shown  for  two  cases:  upper  layer 

stronger  than  the  lower  layer  and  upper  layer  weaker  than  the  lower 
layer.  An  approximate  procedure,  based  on  these  composite  slip 
line  fields,  is  given  for  the  estimation  of  bearing  stresses  in  two- 
layer  soils.  This  approximate  procedure  is  applied  in  a tire-soil 
model  expanded  for  the  consideration  of  two-layer  soils.  The  simu- 
lation of  tire  performance  by  this  expanded  model  is  compared  with 
results  of  small-scale  mobility  tests  performed  at  the  U.S.  Army 
Engineer's  Waterways  Experiment  Station  (WES)  in  layered  soils. 

Field  methods  of  determining  soil  properties  in  layered  soils 
are  examined,  and  a modification  of  ring  shear  tests  is  recommended 
for  the  determination  of  the  strength  of  individual  layers  in 
layered  soils. 
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I.  SCOPE  OF  WORK 


The  scope  of  work,  as  originally  described  in  the  RFP  work 
statement,  was  to  establish  a mathematical  simulation  of  tire- 
soil  interaction  for  soil  conditions  that  vary  with  depth.  With- 
in this  general  scope  emphasis  was  placed  on  the  following  items: 

® Development  of  solution  methods  for  the 
differential  equations  of  plasticity 
for  soil  strength  varying  continuously 
with  depth 

0 Development  of  solution  methods  for  the 
differential  equations  of  plasticity  for 
discrete  variation  of  soil  strength  with 
depth  (layered  systems) 

• Application  of  the  above  solution  methods 
in  tire-soil  interaction  simulation 

® Evaluation  of  present  field  techniques  of 

soil  property  determination  to  assess  their 
suitability  to  characterize  soil  with  a 
variable  strength  profile 
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II.  THEORETICAL  BACKGROUND 


A characteristic  feature  of  soil  deposits  is  that  their 
strength  and  other  properties  inherently  vary  with  depth.  In 
soil  mechanics,  many  theories  and  analytical  methods  have  been 
put  forward  to  allow  proper  consideration  of  this  feature  in  the 
solution  of  the  various  soil  engineering  problems.  As  early  as 
in  1934  Froehlich  (Ref.  1)  proposed  the  use  of  an  elastic  modulus 
linearly  increasing  with  depth  for  the  computation  of  settlements 
and  developed  appropriate  formulas  using  the  theory  of  elasticity. 
Burmister  applied  the  theory  of  elasticity  for  the  computation  of 
stresses  and  displacements  in  layered  systems  (Ref.  2).  In  elas- 
ticity, he  considered  two  kinds  of  boundary  conditions  at  the  layer 
interfaces:  1)  full  continuity  of  stresses  and  displacements 

across  the  interface  (no  limitations  on  shear  stresses  on  account 
of  shear  strength  of  soils)  and  2)  continuity  of  normal  stress 
and  normal  displacement  with  a frictionless  interface  that  allows 
relative  displacement  of  the  adjoining  layers  along  the  interface. 
In  his  theory  each  layer  is  considered  homogeneous  and  is  charac- 
terized by  its  elastic  modulus  and  Poisson's  ratio.  The  numerical 
evaluation  of  the  general  equations  for  layered  systems  set  forth 
by  Burmister  indicated  that  the  boundary  conditions  at  the  layer 
interface,  the  ratio  of  the  elastic  moduli  of  the  adjoining  layers, 
and  the  ratio  of  the  upper  layer  thickness  to  the  dimensions  of  the 
loaded  area  strongly  influence  the  distribution  of  stresses  due  to 
surface  loading. 

In  tire-soil  interaction  problems  of  interest  the  stress 
levels  in  soil  far  exceed  those  for  which  the  theory  of  elasticity 
applies.  In  the  theory  of  plasticity  there  is  no  general  treat- 
ment of  layered  systems.  Sokolovskii  (Ref.  3)  treated  a very 
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special  case  of  layered  systems  — that  of  a "lamellar  medium"  — 
where  the  layering  is  such  that  on  horizontal  planes  the  shear 
stresses  may  not  exceed  a value  defined  by  a Coulomb  type  failure 
condition  along  the  horizontal  plane.  In  the  theory  of  plasticity 
for  soils,  as  it  was  pointed  out  in  Ref.  4,  the  deformations  in 
the  soil  that  occur  prior  to  the  plastic  state  are  disregarded, 
and  the  stress  states  are  calculated  solely  on  the  basis  of  Mohr- 
Coulomb  yield  criterion.  As  a consequence,  in  a logical  expansion 
of  the  plasticity  theory  to  layered  systems  the  boundary  conditions 
at  layer  interfaces  have  to  be  defined  in  terms  of  stresses  alone, 
as  follows:  Normal  and  shear  stresses  across  the  interface  boundary 

must  be  continuous.  The  shear  stress  at  the  interface  may  not  ex- 
ceed the  shear  strength  of  either  adjoining  layer.  In  mathematical 
terms : 


Ti^r  °n  tan  tpi 


t.  < c + o tan  cp 
1 - u n u 


(1) 


In  plasticity  theory,  elastic  deformations  are  not  considered. 
However,  the  displacements  due  to  plastic  flow  are  subject  to  cer- 
tain restraints.  (A  detailed  discussion  of  the  theory  of  velocity 
fields  that  applies  to  these  restraints  is  given  in  Ref.  5.)  The 
application  of  the  theory  of  velocity  fields  to  the  boundary  con- 
ditions at  the  interface  yields  the  following  boundary  conditions: 
Velocity  vectors  normal  to  the  interface  must  either  be  equal  in 
both  layers  or  greater  in  the  upper  layer  than  in  the  lower  layer 
so  that  separation  along  the  boundary  cannot  occur.  Velocity  vec- 
tors tangential  to  the  interface,  however,  can  and  usually  will  be 
different  so  that  one  layer  can  slide  past  the  other. 

Solutions  of  the  plasticity  theory  for  soils  constitute  zones 
of  plastic  equilibrium,  also  known  as  slip  line  fields.  Certain 
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problems  such  as  the  bearing  capacity  or  tire-soil  interaction 
problems  require  a composition  of  various  types  of  zones  (active, 
passive,  and  radial  zones).  Solutions  in  layered  systems  will 
also  be  composed  of  various  zones  that  meet  the  boundary  condi- 
tions at  the  layer  interface.  The  adjoining  zones  in  the  two 
layers  are  not  required  to  be  of  the  same  type,  and  various  con- 
figurations are  possible  as  long  as  the  boundary  conditions  at 
the  interface  are  satisfied.  These  are  discussed  in  more  detail 
in  Section  IV. 

The  consideration  of  a continuous  variation  of  strength  proper- 
ties in  the  solution  of  the  differential  equations  of  plasticity  is 
discussed  in  Ref.  6 for  the  case  of  a nonlinear  Mohr-Coulomb  strength 
envelope.  The  case  where  the  strength  properties  vary  with  depth 
poses  essentially  the  same  type  of  problem  as  a nonlinear  failure 
envelope  and  requires  only  a modification  of  the  numerical  solution 
procedures . 

The  problem  of  bearing  capacity  of  soils  exhibiting  strength 
properties  that  vary  with  depth  either  discretely  or  continuously 
has  been  treated  by  various  researchers,  most  of  them  using  the 
single  failure  surface  mechanism  approach  applied  in  soil  mechanics 
for  stability  analyses.  Tcheng  (Ref.  8)  performed  bearing  tests 
on  layered  systems  using  sand  and  a special  grease  to  model  co- 
hesive soils.  Brown  and  Meyerhof  (Ref.  9)  investigated  the  bearing 
capacity  of  two-layer  purely  cohesive  soils.  Salencon  (Ref.  10) 
treated  the  problem  of  bearing  capacity  of  a purely  cohesive  soil 
with  linearly  varying  strength  analytically.  Yamaguchi  and  Terashi 
(Ref.  11)  presented  analytical  solutions  to  special  bearing  capacity 
problems  on  multilayered  ground. 
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III.  CONSIDERATION  OF  CONTINUOUS  VARIATION  OF  SOIL 


STRENGTH  WITH  DEPTH  IN  THE  TIRE -SOIL  MODEL 


The  methods  of  numerical  computations  for  the  determination 
of  the  geometry  of  slip  line  fields  and  associated  stresses  have 
been  discussed  in  detail  in  Refs.  4 through  6.  For  brevity,  only 
the  differences  in  the  computation  methods  necessitated  by  the 
consideration  of  strength  variation  with  depth  are  discussed  herein. 

The  variation  of  the  strength  properties  with  depth  considered 
in  the  tire-soil  model  are  shown  in  Figs.  1 and  2.  Figure  1 shows 
the  case  where  the  f intercept  at  the  cr-axis  is  constant,  while 
Fig.  2 shows  the  more  general  stiuation  when  the  f intercept 
varies  with  the  assumed  variation  of  soil  strength.  Because  of 
this  variation  of  the  f intercept  the  numerical  solution  methods 
for  the  differential  equations  of  plasticity  for  soils  presented 
in  Ref.  6 have  to  be  modified.  In  these  solutions  the  principal 
stress  variable,  a,  is  defined  as  the  distance  of  the  center  of 
the  Mohr  circle  from  the  intersection  of  Mohr-Coulomb  envelope  with 
the  a~axis  (-^) . If  rjj  varies  with  depth,  then  the  o values 
have  to  be  adjusted  to  this  variation.  Figure  3 shows  the  schematic 
diagram  of  the  computation  of  variables  at  an  i,j  nodal  point  with 
respect  to  this  variation.  In  practice,  it  was  found  that  one 
iteration  on  the  f values  for  the  determination  of  the  appro- 
priate strength  values  was  satisfactory. 

A more  detailed  flow  chart  for  the  computation  of  a single 
slip  line  field  (subroutine  SLFI)  for  variable  strength  is  given 
in  the  appendix.  In  the  computer  program  represented  by  the  flow 
chart,  the  following  variation  of  the  strength  parameters  with 
depth  is  assumed. 
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Fig.  1 Strength  Variation  with  Depth  z for  Constant  ip  Intercept 


Fig.  2 Strength  Variation  with  Depth  z Without  Restrictions 
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Fig.  3 Scheme  for  Computing  the  Variables  a and  at  Nodal 
Point  i,  j for  Strength  Varying  with  Depth  z. 
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Designation  of  the  variables  in  the  subroutine  for  the  com- 
putation of  a single  slip  line  field  is  as  follows.  Double  sub- 
scripted variables: 


X(I,J)  = x 

X J J 

Z(I,J)  = z 

5 J 

S(I,J)  = a.  . 

J J 

T(I , J)  = 0.  . 
P(I,J)  - a.  . 


(x  coordinate  of  nodal  point) 

(z  coordinate  of  nodal  point) 

( i(<^  + <73)  + f) 

(direction  of  major  principal  stress) 

^ °l  + a3> 


The  dimension  statements  for  these  subscripted  variables  cor- 
respond to  a 43  x 15  grid  (I  = 1 to  45,  J = 1 to  15)  for  the  com- 
putation of  the  geometry  of  slip  line  fields. 


Single  subscripted  variables  designate  the  following  values  of 
the  variables  at  the  interface  at  the  "jn  line  intersections. 


D(J)  - 
H(J)  = 
Q(J)  - 
E(J)  - 
U(J)  = 
V(J)  = 
A(J) , B(J), 


a . 
J 

qj 


r . 
J 

x. 

J 


(interface  friction  angle) 
(central  angle) 

(normal  stress) 

(shear  stress) 

(x  coordinate  of  interface) 
(z  coordinate  of  interface) 


C(J)  are  auxiliary  variables. 
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Constants  and  Parameters: 


R0 

= radius 

G0 

= y 

(unit  weight) 

COH 

- C 

(cohesion  at  surface) 

PHI 

= <P 

(friction  angle  at  surface) 

Dl 

= 5 

(interface  friction  angle) 

A0 

= a 

e 

(entry  angle  in  arc) 

AM 

= a 

m 

(central  angle  at  end  of  field) 

Al 

= a 

s 

(tangent  angle  to  log  spiral) 

PL 

= Pi 

(limit  pressure) 

XX 

= - 1 

for  front  field 

XX 

= + 1 

for  rear  field 

XX 

= 2 

for  rear  field  only 

DEPHI 

■ \ 

(constant  for  9 variation  with  depth) 

DEC  ON 

= Dcoh 

(constant  for  c variation  with  depth) 

The  main  program  for  the  computation  of  drawbar  pull,  torque, 
and  sinkage  is  essentially  the  same  as  reported  for  the  tire-soil 
model  for  uniform  strength  in  Ref.  7. 

Typical  results  obtained  with  the  tire-soil  model  modified 
for  variable  strength  are  shown  in  Fig.  4 where  pull  performance 
relationships  are  shown  for  four  cases,  each  depicting  a different 
variation  of  soil  strength  with  depth.  Common  input  conditions 
for  these  cases  were: 
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PULL  COEFFICIENT 


SLIP  % 

Fig.  4 Comparison  of  Pull  Performances  for  Soil  Strength 
Varying  with  Depth 
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Tire  Characteristics 


Tire  size 

9.00 

- 14 

Nominal  radius 

1.18 

ft 

Nominal  width 

0.74 

ft 

Inflation  pressure 

18.7 

psi 

Deflection 

25% 

Deflection  coefficient 

0.88 

Load 

850  : 

lb 

Yuma 

Sand 

Case  No.  1 represents  uniform  conditions  corresponding  to  a 
cone  index  gradient  (CGR)  of  15,  or  cp  = 44.7°.  In  Case  No.  2, 
the  soil  strength  decreases  linearly  from  cp  = 44.7°  at  the  sur- 
face to  cp  = 39.7°  at  the  deepest  point  of  the  slip  line  field. 

Case  No.  3 again  represents  uniform  conditions  where  the  soil 
strength  equals  cp  = 39.7°,  the  minimum  soil  strength  in  Case 
No.  2.  This  friction  angle  corresponds  to  a CGR  of  5.5.  Finally, 
in  Case  No.  4,  the  soil  strength  increases  linearly  from  cp  = 39.7° 
to  cp  = 44.7°  at  the  deepest  point  of  the  s.lip  line  fields.  As 
is  seen  from  Fig.  4,  the  pull  performance  for  variable  shear  strength 
is  close  to  that  in  a uniform  soil  having  the  same  strength  as  the 
variable  strength  profile  at  the  surface. 

Figure  5 shows  the  computed  distribution  of  normal  and  shear 
stresses  and,  in  the  upper  part,  the  deflection  of  the  tire  for 
Case  No.  1 at  20  percent  slip,  as  obtained  with  the  computer 
graphics  program.  Figure  6 shows  the  same  for  Case  No.  2.  Note 
that  in  the  radial  plot  a rigid  wheel  with  the  nominal  tire  diame- 
ter is  shown  for  the  convenience  of  plotting;  the  tire  deformation 
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TIRE 

DEFORMATION 


o = NORMAL  STRESS  IN  LB/SQ  FT 
T = SHEAR  STRESS  IN  LB/SQ  FT 


Fig.  5 Distribution  of  Normal  and  Shear  Stresses  and  Tire 
Deflection,  Case  No.  1,  20 % Slip 
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Fig. 


TIRE 

DEFORMATION 


a = NORMAL  STRESS  IN  LB/SQ  FT 
t = SHEAR  STRESS  IN  LB/SQ  FT 


6 Distribution  of  Normal  and  Shear  Stresses  and  Tire 
Deflection,  Case  No.  2,  20%  Slip 
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is  shown  separately  at  the  top.  The  difference  in  the  interface 
stresses  computed  for  the  two  cases  is,  indeed,  minor  and  results 
in  little  difference  in  pull  performance. 

The  effect  of  a decreasing  shear  strength  on  the  slip  line 
fields  that  develop  underneath  a tire  is  that  the  slip  line  fields 
encompass  a greater  depth  than  those  for  uniform  strength. 

Figure  7a  shows  a typical  rear  slip  line  field  obtained  for 
the  conditions  of  Case  2 where  the  shear  strength  decreased  from 
about  45°  at  the  surface  to  about  33°  at  the  point  of  the 
field.  For  comparison,  the  slip  line  field  for  Case  1 uniform 
strength  (9  = 45°)  is  also  shown  (Fig.  7b). 

The  growth  of  the  slip  line  field  in  the  z direction  due  to 
a decreasing  shear  strength  results  in  encountering  layers  of  less 
and  less  strength.  At  some  rate  of  strength  decrease,  a progressive 
process  takes  place  for  which  no  slip  line  field  solution  may  be 
found.  This  is  interpreted  as  a condition  where  complete  and  deep 
reaching  failure  would  take  place  in  the  soil  under  the  tire  load, 
indicating  a “no  go"  condition. 

In  conclusion,  in  sand  the  surface  strength  is  the  dominating 
factor,  but  there  is  a limit  to  the  rate  of  decrease  at  which  a 
no  go  situation  occurs.  In  cohesive  soils  the  situation  is  very 
similar,  as  far  as  the  limit  of  strength  decrease  is  concerned  at 
which  no  go  condition  sets  in.  The  effect  of  strength  decrease  on 
pull  performance  is  less  pronounced  than  in  sand  and  sometimes  para- 
doxical in  that  higher  pull  performance  is  obtained  with  decreasing 
strength.  The  explanation  of  this  paradox  is  shown  in  Fig.  8.  The 
development  of  the  interface  shear  stress  depends  on  the  angle  b 
and  equals 


b)  FOR  CASE  1,  4>  = 45° 


Fig.  7 Rear  Slip  Line  Fields 
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In  this  equation  tan  6 is  the  contribution  of  the  cohesion 
to  the  interface  stress  that  may  be  called  adhesion.  The  pull  per- 
formance is  governed  by  the  ratio  of  the  interface  shear  stress  to 
the  normal  stress,  or  tan  6'  in  Fig.  8.  A decreasing  shear 
strength  reduces  the  normal  stresses  while  not  affecting  t that 
depends  on  the  surface  shear  strength  and  6.  Thus,  in  cohesive 
soils  the  T/an  ratio  generally  increases  as  the  shear  strength 
decreases  with  depth.  If  the  effect  of  the  increase  of  the  T/a 

n 

ratio  is  not  set  off  by  an  increase  in  sinkage,  then  pull  perfor- 
mance improves  with  a shear  strength  that  decreases  with  depth. 


Fig.  8 Development  of  Interface  Shear  Stresses  in  Cohesive  Soil 

The  conclusions  that  can  be  drawn  from  the  tire-soil  model  for 
soil  strength  varying  with  depth  is  that  the  surface  strength  has 
more  influence  on  the  pull  performance  than  the  strength  of  the 
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underlying  layers  as  long  as  the  strength  decrease  does  not  result 
in  a critical  no  go  situation.  It  is  recommended  that  criteria 
for  the  rate  of  strength  decrease  that  results  in  such  a situation 
be  established  by  further  research. 

In  these  analyses  a linear  variation  of  soil  strength  with 
depth  was  assumed.  Other  relationships  that  describe  strength 
variation  with  depth  may  also  be  assumed  and  appropriately  intro- 
duced into  the  program.  No  difficulties  are  expected  to  occur 
with  other  forms  of  strength  variation  as  long  as  the  variation 
is  smooth. 


17 


IV.  THE  EFFECT  OF  DISCRETE  LAYERING  ON  THE  MAGNITUDE 
AND  DISTRIBUTION  OF  FAILURE  STRESSES 


The  concept  of  tire-soil  interaction  developed  for  homogene- 
ous soil  conditions  postulates  that  the  soil  is  in  the  plastic 
state  of  failure  whenever  the  failure  stresses  at  the  tire-soil 
interface  are  less  than  the  limit  pressure  characteristic  of  the 
tire.  Stresses  in  soil  in  the  plastic  state  of  failure  are  de- 
termined by  the  numerical  integration  of  the  differential  equa- 
tions of  plasticity  for  soils;  the  orientation  of  failure  surfaces 
and  the  extent  of  the  plastic  state  are  represented  by  the  slip 
line  field  obtained  in  the  course  of  numerical  integration. 

To  develop  a method  for  determining  the  effect  of  discrete 
layering  on  tire-soil  interaction,  the  problem  of  constructing  a 
single  slip  line  field  must  first  be  solved.  Since  this  problem 
has  not  been  treated  previously  in  trwo-layer  soils,  it  is  expedient 
to  choose  an  approach  where  the  problems  of  discrete  layering  may 
be  solved  without  the  complications  inherent  in  the  slip  line  field 
calculations  in  the  tire-soil  interaction  problem.  The  simple  two 
dimensional  bearing  capacity  case  is  suitable  for  this  purpose. 

PROBLEMS  ENCOUNTERED 

Originally  it  was  planned  to  develop  solution  methods  for  the 
bearing  capacity  of  two-layer  soils  by  using  a common  grid  system 
in  both  layers  and  use  auxiliary  nodal  points  at  the  locations 
where  the  slip  lines  cross  the  interface  boundary.  Computer  pro- 
grams were  written  for  the  necessary  algorithms  for  the  computation 
of  the  coordinates  and  stress  states  at  these  auxiliary  points,  but 
unexpected  problems  arose  in  connection  with  certain  situations 
illustrated  in  Fig.  9.  In  the  numerical  computation  scheme,  the 
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Fig.  9 Computation  of  Nodal  Points  at  Layer  Interfaces 
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variables  are  computed  at  an  i,j  point  from  the  values  at  the 
i-l,j  and  i,j-l  points.  If  both  of  these  points  are  in  the 
same  layer  (as  shown  in  Fig.  9a),  there  is  no  problem.  If  the 
two  points  are  in  different  layers,  but  hypothetical  points  (as 
they  were  in  the  adjacent  layer)  may  be  computed  (as  shown  in 
Fig.  9b),  the  problem  solution  requires  additional  storage  of  data 
for  hypothetical  points.  The  problem  for  which  no  satisfactory 
solution  has  been  found  is  shown  in  Fig.  9c, where  the  nodal  points 
required  to  compute  a hypothetical  point  for  i,j  are  not  in  the 
same  layer  and  a hypothetical  i,j-l  point  is  not  available  from 
previous  computations.  Various  methods  have  been  tried  to  over- 
come this  problem  by  using  estimated  values  for  these  points,  but 
none  has  proved  to  be  generally  satisfactory. 

From  the  many  trial  computations  performed  for  this  approach 
it  became  evident  that  using  separate  grid  systems  in  the  two  layers 
has  a better  chance  of  success.  It  was  also  concluded  that  dif- 
ferent solution  methods  would  be  required,  depending  upon  whether 
the  upper  layer  is  stronger  than  the  lower  one,  or  vice  versa. 

RELATIVE  STRENGTH  OF  LAYERS 

There  are  several  cases  with  respect  to  the  relative  strength 
of  the  layers  in  a two-layer  system. 

Case  1 

For  a given  normal  stress,  an>  and  for  cr^  < i.e., 

c + a tan  cp  < c,  + o tan  cp,  (4) 

un  u ln  1 

The  subscripts  u and  1 refer  to  upper  and  lower  layer,  respec- 
tively. This  condition  is  represented  schematically  by  the  Mohr- 
Coulomb  construction  shown  in  Fig.  10. 
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Case  2 


For  a given  normal  stress,  crn,  and  for  < au,  i.e., 

c,  + o tan  cp.  < c + a tan  cp  (5) 

In  lun  u 

This  condition  is  represented  schematically  in  Fig.  10a  also,  ex- 
cept that  the  u and  1 subscripts  are  everywhere  interchanged. 

Case  3 

There  is  a •'crossover"  normal  stress,  a , where  the 

c 

strength  of  the  upper  layer  equals  that  of  the  lower  layer.  For 

o < a , 

n c 

c,  + o tan  cp.  < c + o tan  cp  (6) 

In  l un  u v ' 

and  for  a given  normal  stress,  a , greater  than  the  crossover 

normal  stress,  o , 

c 

c + a tan  cp  < c,  + a tan  cp.  (7) 

un  u ln  1 v 

This  condition  is  shown  schematically  in  Fig.  11. 

Case  4 

Case  4 is  the  reverse  of  Case  3 with  the  u and  1 sub- 
scripts interchanged. 

If  the  maximum  normal  stress  at  the  interface  is  less  than 
ac  in  a certain  problem,  then  Cases  3 and  4 revert  to  Cases  1 and 
2,  respectively.  If  a is  within  the  range  of  normal  stresses 
that  are  transmitted  through  the  interface  in  a certain  situation, 
then  there  is  relatively  little  difference  in  the  strength  of  the 
two  layers,  and  an  acceptable  approximation  may  be  found  by  treat- 
ing the  case  at  homogeneous  soil  with  strength  properties  averaged 
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between  the  two  layers.  For  these  reasons  the  solution  method 
discussed  in  the  following  are  restricted  to  Case  1 and  2. 


SOLUTION  METHOD  FOR  CASE  1 

A solution  of  the  simple  two  dimensional  bearing  capacity 
problem  must  satisfy  the  boundary  conditions  at  both  the  surface 
and  at  the  layer  interface.  The  boundary  conditions  at  the  layer 
interface  require  the  continuity  of  normal  and  shear  stresses 
across  the  interface  with  the  additional  condition  that  the  shear 
stress  may  not  exceed  the  strength  of  the  weaker  layer.  Since  the 
strength  of  soils  depends  on  the  normal  stress,  it  is  convenient 
to  express  the  condition  on  shear  strength  in  terms  of  the  layer 
interface  friction  angle  8 (a  concept  similar  to  the  tire-soil 
interface  friction  angle) . The  strength  of  the  upper  layer  may  be 
expressed  as 

s = (c  + a ) tan  cp 
u v u n 

and  the  shear  stress  at  the  interface 


T 


U 


+ a ) 

n' 


tan  6 

u 


Ti  ■ (ci  + °n>  tan  V 

(cu  + <*r)  tan  6U  = (c^  + ffn)  tan  5 £ < (cu  + an)  tan  cpy 

6 < cp 

u 

c + a c + a 

tan  6 . = ~ - . — ^ tan  6 < ■ — - tan  cp  (8) 

i c , + a u-c.+ct  u v/ 

in  in 


The  shear  stresses  and  6y  angles  generated  at  the  layer  in- 
terface by  a slip  line  field  in  the  upper  layer  vary  from  0 in 
the  passive  zone  to 
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and 


T 

max 


(c  + a ) tan  cp 
v u n'  u 


at  the  point  where  the  interface  is  tangent  to  a "j"  slip  line. 

At  this  point,  the  interface  is  a slip  line  and  = cpu«  The 

variation  of  6y  angle  in  a slip  line  field  is  shown  in  Fig.  12. 

The  slip  line  field  in  the  area  (ABODE)  may  be  constructed  without 

interference  from  the  lower  layer.  The  dash-dot  line  to  the  right 

from  the  point  where  6 =6  indicates  the  variation  of  6 for 

u max 

hypothetical  uniform  soil  conditions,  i.e.,  that  would  occur  along 
the  line  of  the  interface  if  soil  conditions  did  not  change  across 
that  line. 

The  d angle,  the  direction  of  the  major  principal  stress, 
is  computed  from  the  6 angle  as  follows  (Ref.  12) 


(1  - k)  f + (kA  - 6) 


(9) 


where 


A = arc  sin 


sin  6 


sm  cp 

k = +1  in  the  passive  case 
k = -1  in  the  active  case 

The  passive  case  k = +1  includes  all  stress  states  where 

a r > a > a, 
nf  n 3 

and  the  active  case  (k  = -1)  all  stress  states  where 


°nf  < an  < al 


j f = normal  stress  at  the  failure  plane  (Fig.  13) 
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Point  D at  the  interface  (Fig.  12)  (where  o = crn^  and 

au  = :i'S  t*ie  Po:i-nt:  that  separates  the  k = +1  and  k = -1 

stress  states.  The  variation  of  the  6 angle  along  the  line  of 

the  interface  in  the  hypothetical  uniform  case  would  be  continuous 

even  though  k could  change  at  point  "d"  from  -1  to  +1,  since 

Eq.  (9)  yields  the  same  value  for  6 with  either  value,  if  5u  = 

cp  . In  this  case 
u 

A = arc  sin  1 = \ ir 


k = +1  , 

e = -|(^  tt  - cp)  = 

£ TT  - i Cp 

(10) 

k = -1  , 

6 = \ r + i ir 

-<p)“£tt-£<p 

In  the  case  of  the  bearing  capacity  problem,  the  only  con- 
ceivable solutions  are  those  where  the  character  of  stress  state 
(the  "k"  value)  does  not  change  across  the  interface.  (A  change 
in  the  k value  would  result  in  a permissible,  but  inappropriate 
stress  discontinuity  along  the  interface.)  Thus,  the  boundary  con- 
ditions at  the  interface  for  a slip  line  field  solution  in  the  lower 
layer  are  characterized  by  a k = -1  along  (DC)  and  a k = +1 
value  along  (DF) . The  0 values  along  the  interface  can  be  com- 
puted by  Eq.  (9).  At  an  infinitesimal  distance  to  the  right  of 
point  D,  0 in  the  lower  layer  could  be  computed  using  k = +1, 
while  at  an  infinitesimal  distance  to  the  left,  it  could  be  com- 
puted using  k = -1.  While  in  the  upper  layer  this  change  of  the 
sign  of  k does  not  affect  the  continuity  of  0,  in  the  lower 
layer,  where  6 < cp,  there  will  be  a discontinuity  in  0 at  this 
point.  This  situation  has  been  recognized  as  identical  to  that  at 
the  singular  point  of  the  bearing  capacity  case  (point  A in  Fig.  12) 
requiring  a radial  zone  to  accommodate  this  discontinuity.  Since 
the  boundary  conditions  along  (DC)  are  also  defined  by  the  con- 
tinuity condition  on  the  normal  and  shear  stresses,  a slip  line 
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field  can  be  constructed  in  the  lower  layer  for  any  assumed  varia- 
tion of  6 (or  0)  along  DF.  In  lieu  of  anything  better,  it  was 
assumed  that  6 along  DF  would  change  the  same  way  as  in  the  case 
of  uniform  soil  conditions. 

The  slip  line  field  in  the  lower  layer  yields  the  normal  and 
shear  stresses  at  the  DF  portion  of  the  interface  from  which  the 
values  of  all  four  variables  (x,  z,  cr,  and  6)  for  the  upper 
layer  strength  parameters  can  be  computed.  These  constitute  the 
boundary  conditions  for  which  an  active  zone  type  slip  line  field 
may  be  constituted  in  the  upper  layer.  Thus,  the  solution  in  Case  1 
consists  of  three  fields,  as  shown  in  Fig.  14.  Field  (abed)  is  a 
regular  slip  line  field  cut  off  at  the  interface.  Field  (edefghi) 
is  a bearing  capacity  field  consisting  of  the  regular  three  zones 
but  for  a variable  6 as  boundary  conditions  at  the  interface. 

Field  (efij)  consists  only  of  an  active  zone.  For  the  computation 
of  this  field,  auxiliary  points  are  needed  in  the  zone  bounded  by 
the  line  (jklm). 

A computer  program  was  written  for  the  computation  of  a com- 
posite slip  line  field  comprising  the  three  types  of  fields  shown 
in  Fig.  14.  Typical  composite  slip  line  fields  are  shown  in  Fig.  15a 
through  15g  for  the  following  strength  properties  of  the  two  layers: 


Upper  Layer  | 

Lower  Layer 

■ 

USM 

EM 

RHI 

■H 

3HI 

SKSSIHKIui 

BBBSm 

tSSRSft 

15a 

10 

30 

100 

14.50 

40 

100 

15b 

10 

30 

100 

10 

35 

100 

15c 

80 

20 

100 

100 

25 

100 

15d 

50 

20 

100 

100 

20 

100 

15e 

50 

10 

100 

80 

10 

100 

15f 

10 

20 

100 

10 

30 

100 

15g 

10 

10 

100 

20 

12 

100 
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The  construction  of  these  composite  slip  line  fields  is  not 
necessarily  feasible  for  all  combinations  of  the  strength  proper- 
ties of  the  two  layers.  If  the  upper  layer  is  very  weak  compared 
to  the  lower  one,  then  the  last  portion  of  the  composite  field 
(eijk  in  Fig.  14)  cannot  be  constructed.  The  reason  for  this  is 
that  the  upper  weak  layer  cannot  transmit  stresses  high  enough  to 
result  in  a plastic  state  of  failure  in  the  lower,  stronger  layer. 
The  composite  field  represents  situations  where  both  layers  are  in 
the  plastic  state  of  failure  under  the  applied  load. 

The  case  when  the  relative  strength  of  the  two  layers  is  such 
that  the  lower,  stronger  layer  is  not  in  the  plastic  state  failure 
under  the  applied  load  requires  another  method  of  analysis  of  the 
stresses  in  the  upper  layer.  Generally,  the  constraints  on  a stress 
field  in  the  upper  layer  exercised  by  a rigid  base  are  such  that 
there  is  no  complete  solution  for  the  bearing  capacity  problem  (an 
incomplete  solution  would  leave  the  stress  states  in  certain  areas 
unresolved).  In  such  a case  kinematically  admissible  slip  line 
fields  (upper  bound  solutions)  may  be  used  to  estimate  an  upper 
bound  for  the  bearing  stresses.  Such  a field  for  the  case  of  cp  = 0 
is  shown  in  Fig.  16.  Implications  of  such  a slip  line  field  are 
discussed  in  connection  with  the  problem  of  slipperiness  in  Sec- 
tion VI. 

Typical  bearing  stresses  calculated  on  the  basis  of  composite 
slip  line  fields  are  shown  in  Figs.  17a  through  17e  for  the  fol- 
lowing conditions: 
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LAYER  2 


Fig.  16  Kinematic  Slip  Line  Field  for  <j>  = 0,  Strong 
Bottom  Layer 


In  these  figures  square  or  diamond  symbols  indicate  the 
bearing  stresses  obtained  from  the  composite  slip  line  fields. 
Dashed  lines  show,  for  comparison,  bearing  stresses  computed  for 
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17  Bearing  Stresses  Computed  on  the  Basis  of  Composite 
Slip  Line  Field  for  Two-Layer  Systems,  Sheet  1 of  2 
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Fig.  17  Bearing  Stresses  Computed  on  the  Basis  of  Composite 
Slip  Line  Fields  for  Two-Layer  Systems,  Sheet  2 of  2 


hypothetical  uniform  soil  conditions,  assuming  that  the  uniform 
soil  has  the  strength  properties  of  either  the  upper  or  the  lower 
layer. 

Figures  17a  through  17e  indicate  that  there  is  a jump  in  the 
normal  stresses  at  the  point  where  the  lower  layer  starts  to  af- 
fect the  slip  line  field  in  the  upper  layer.  It  is  necessary  to 
examine  whether  such  a jump  is  permissible  or  not . One  way  to 
analyze  this  situation  is  to  consider  the  normal  stresses  to  the 
left  of  the  jump  as  a surcharge  and  analyze  whether  a bearing 
failure  due  to  the  jump  in  the  normal  stresses  is  possible.  In 
the  cases  shown  the  stress  difference  at  the  jump  is  not  enough  to 
cause  such  a failure,  indicating  that  the  upper  layer  is  able  to 
sustain  the  stress  difference.  Nevertheless,  the  slip  line  field 
associated  with  such  stress  differences  is  a different  one  from 
that  of  the  composite  field  and  one  could  speculate  that  the  stress 
states  in  the  soil  would  adjust  to  minimize  the  stress  jump.  Such 
an  adjustment  may  take  place  by  the  development  of  opposing  shear 
stresses  beneath  the  footing.  These  shear  stresses  could  increase 
the  lower  bearing  stresses  to  the  right  of  the  jump  and  decrease 
the  higher  bearing  stresses  to.  the  left.  The  effect  of  such  shear 
stresses  is  illustrated  in  Fig.  17f,  which  shows  the  bearing  stresses 
computed  for  an  interface  friction  angle  of  6 = +5°  from  the  edge 
of  the  footing  to  the  stress  jump  and  for  8 = -5°  past  the  jump. 

The  figure  shows  that  these  opposite  shear  stresses  tend  to  lessen 
the  jump  of  the  magnitude  in  the  normal  stresses  and  could  con- 
ceivably completely  eliminate  the  jump.  The  composite  slip  line 
field  for  these  conditions  is  shown  in  Fig.  18. 

This  equalization  of  the  normal  stresses  due  to  the  develop- 
ment of  shear  stresses  is  very  likely  to  occur  in  actuality.  Un- 
fortunately, there  is  neither  a theoretical  basis  nor  experimental 
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information  that  would  allow  to  estimate  these  shear  stresses  or 
their  distribution  quantitatively.  Because  of  this  unknown  factor, 
the  bearing  stresses  computed  without  the  consideration  of  shear 
stress  developments  are  considered  as  limits  only.  Further  theo- 
retical and  experL mental  research,  including  measurements  of  the 
shear  stresses  under  various  conditions,  is  needed  for  better  under- 
standing of  the  phenomena  and  the  development  of  a theory  to  take 
these  shear  stresses  properly  into  account. 


SOLUTION  METHOD  FOR  CASE  2 


In  Case  2 the  upper  layer  is  the  stronger  one.  Consequently, 
the  strength  of  the  lower  layer  limits  the  shear  stresses  that  are 
generated  in  an  upper  layer  slip  line  field  along  the  interface. 
The  construction  of  an  upper  layer  slip  line  field  cannot  be  con- 
tinued if  at  any  location  along  the  interface  the  following  condi- 
tion holds: 


T > T 


max 


(cf  + an)  tan  ^ 


(11) 


If  the  limitation  on  the  shear  stresses  is  ignored  and  the 
construction  of  the  slip  line  field  continued,  hypothetical  shear 
stresses  along  the  interface  may  be  determined.  These  indicate 
that  the  shear  stresses  exceed  the  shear  strength  of  the  lower 
layer  along  some  length  of  the  interface.  Thus,  this  length  of  the 
interface  must  be  a slip  line  in  a field  in  the  lower  layer,  a basic 
feature  of  a failure  mechanism  that  is  intuitively  recognized  as 
correct.  This  feature  has  to  be  incorporated  in  a composite  slip 
line  field  that  comprises  the  plastic  zones  in  both  layers.  Fig- 
ure 19  shows  schematically  the  various  zones  of  such  a composite 
slip  line  field. 
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Fig.  18  Composite  Slip  Line  Field  With  Base  Friction 


Fig.  19  Scheme  for  Composition  of  Slip  Line  Fields  for 
Case  2 


In  the  upper  layer,  field  (abed)  is  composed  of  the  usual 
three  zones  (active,  passive,  and  radial).  The  field  is  bounded  at 
the  left  side  by  the  "i"  line,  which  at  the  crossover  point  at 
the  interface  produces  shear  stresses  equal  to  the  shear  strength 
of  the  lower  layer.  In  the  lower  layer,  the  field  (cdefgh)  is  a 
special  one,  which,  contrary  to  the  usual  fields,  does  not  have  a 
radial  zone  and  a singular  point.  Instead,  a "square"  zone  provides 
the  transition  from  the  active  zone  to  the  passive  zone.  The  top 
line  of  this  square  zone  is  a "j"  slip  line.  To  construct  this 
field,  the  6 values  along  (ef)  are  assumed  to  vary  in  a similar 
manner  as  they  would  in  a field  for  uniform  soil  conditions.  The 
field  (efij)  is  similar  to  that  in  Case  1.  Auxiliary  points  are 
needed  (within  the  area  indicated  by  dashed  line)  for  the  computa- 
tion of  this  field. 

The  method  of  computation  of  each  of  these  fields  is  essentially 
the  same  as  in  Case  1 except  for  the  square  field  (degh) , which  re- 
quires special  consideration. 

The  top  j line  in  this  field  coincides  with  the  interface 
where  the  interface  friction  angle  6 equals  cp  and  is  constant 
along  (de) . The  angle  of  inclination  of  the  major  principal  stress 
(cr^)  is  also  constant  and  equals 

e = -i  7T  - i cp  (12) 

For  the  computation  of  this  field,  values  of  all  four  variables 
(x,  z,  ct,  and  6)  of  the  differential  equations  of  plasticity  must 
be  given  along  (de) . At  point  d the  normal  and  shear  stress 
transmitted  from  the  upper  layer  are  known,  and  the  value  of  a 
may  be  computed  therefrom.  The  distance  (de)  is  divided  into 
equal  intervals  designating  the  end  locations  of  i lines.  For 
the  determination  of  the  a values  associated  with  these  points, 
the  following  considerations  apply. 
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The  differential  equation  of  plasticity  that  governs  the  a 
variation  is  as  follows  (Ref.  3): 

dcr  + 2a  tan  cp  d0  = ^(dz  + tan  cp  dx)  (13) 

As  explained  previously,  along  (de)  6 = constant,  and  for 
an  interface  parallel  with  the  x-axis,  z = constant.  Conse- 
quently d0  = 0 and  dz  = 0 and  Eq.  (13)  is  reduced  to 

da  = y tan  <p  dx  (14) 

This  allows  the  computation  of  a along  (de) . 

With  the  values  of  all  four  variables  defined  along  (de) , 
the  slip  line  field  in  the  square  zone  may  be  computed.  The  length 
of  (de)  is,  however,  not  known  and  has  to  be  estimated.  The  com- 
putation of  the  rest  of  the  composite  field  is  similar  to  that  in 
Case  1. 

Figures  20a  and  20b  show  typical  composite  slip  line  fields 
for  the  following  conditions: 
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The  computed  bearing  stresses  for  these  cases  are  shown  in 
Figs.  21a  and  21b.  The  jump  in  the  bearing  stresses  shown  in  these 
figures  would  be  alleviated  by  opposite  directed  shear  stresses  de- 
veloping at  the  base  of  the  footing,  just  as  it  was  shown  for 
Case  1. 
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EVALUATION  OF  THE  EFFECTS  OF  DISCRETE  LAYERING  ON  THE  FAILURE 
NORMAL  STRESSES 


The  distribution  of  bearing  stresses  obtained  from  the  com- 
posite slip  line  fields  for  a two-layer  system  is  characterized  as 
consisting  of  two  parts  separated  by  a jump.  The  first  part  ex- 
tends from  the  edge  to  a distance  within  which  a slip  line  field 
for  the  upper  layer  can  be  constructed  without  interference  from 
the  lower  layer.  The  second  part  extends  from  this  point  to  the 
end  of  the  composite  slip  line  field.  In  the  first  part  the  bear- 
ing stresses  are  identical  with  that  for  uniform  soil  with  the 
strength  of  the  upper  layer.  In  the  second  part  the  bearing  stresses 
closely  parallel  that  for  uniform  soil  with  the  strength  of  the  lower 
layer,  but  there  is  a shift  in  the  magnitude  of  these  stresses.  For 
Case  1 this  shift  in  the  bearing  stresses  can  be  approximately  deter- 
mined, as  shown  in  Fig.  22.  The  bearing  stress  at  the  point  of  in- 
terference ("A")  (where  the  j slip  line  ending  at  this  point 
touches  the  layer  interface)  equals  that  for  the  first  layer,  plus 
the  difference  between  the  edge  bearing  stress  for  upper  layer 
strength  and  that  for  lower  layer  strength  computed  as  if  the  soil 
were  uniform. 


In  Case  2,  the  bearing  stress  at  this  point  can  be  approxi- 
mately calculated  from  the  a value  at  point  "e"  (Fig.  19), 
since  the  a value  is  practically  unchanged  in  the  field  (efij). 
The  normal  stresses  at  (ij)  are,  however,  affected  by  the  6 
value  at  (ij)  and  can  be  calculated  as  follows: 


q = o cos  5 (cos  6 + 


l 


2k 
cos  6 


- cos  cp)  - f 


(15) 


Thus,  the  bearing  stresses  in  a two-layer  system  can  be  ap- 
proximately calculated  from  the  bearing  stresses  calculated  for 
uniform  soil  with  the  strength  of  each  of  the  layers  and  from  the 
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BEARING  STRESSES  CALCULATED 
WITH  LAYER  2 STRENGTH  PARAMETERS 


j LAYER  2 

t 


Fig.  22 


Approximate  Determination  of  Bearing  Stresses 
for  Two-Layer  Soils 


1 


magnitude  of  the  shift  at  point  A calculated  by  the  methods  dis- 
cussed above.  These  bearing  stresses,  however,  can  be  considered 
as  limits  only.  The  actual  bearing  stress  depends  on  the  develop- 
ment of  oppositely  directed  shear  stresses  at  the  base.  These 
shear  stresses  reduce  the  magnitude  of  the  shift  in  the  bearing 
stresses . 
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V.  TIRE-SOIL  MODEL  FOR  TWO-LAYER  SOILS 


The  methods  of  constructing  composite  slip  line  fields  (de- 
scribed in  Section  IV)  for  the  bearing  capacity  problem  in  two- 
layer  soils  can  be  directly  applied  in  the  tire-soil  model.  How- 
ever, the  resulting  discontinuous  stress  distribution  would  present 
serious  problems  in  the  various  iteration  schemes  employed  in  that 
model.  For  example,  the  entry  angle  there  is  found  by  iteration  on 
the  condition  that  the  normal  stress  at  the  angle  of  separation, 
qm,  computed  from  the  slip  line  field,  should  match  the  limit 
pressure  p^  within  an  allowable  tolerance  limit.  If  the  value 
of  qm  falls  within  the  magnitude  of  the  jump  of  the  normal  stresses 
at  the  point  of  stress  discontinuity,  then  this  condition  cannot  be 
satisfied  and  some  other  means  of  finding  the  entry  angle  would  have 
to  be  devised.  At  any  rate,  the  more  complex  computations  associated 
with  the  construction  of  composite  slip  line  fields  and  the  problems 
of  iterations  in  the  case  of  stress  discontinuity  would  require  a 
very  elaborate  computation  scheme  accompanied  by  a great  increase  in 
computing  time. 

The  methods  for  estimating  the  interface  normal  stresses  in  a 
two-layer  system  on  the  basis  of  normal  stresses  computed  for  uni- 
form soil  conditions  (described  in  Section  IV)  offer  a more  practi- 
cal solution  for  the  tire-soil  interaction  problem  in  two-layer 
soils.  Instead  of  using  composite  slip  line  fields  for  the  compu- 
tation of  interface  stresses,  the  approximate  method  is  applied  for 
the  estimation  of  these  stresses.  Accordingly,  the  interface  stresses 
are  computed  first  for  the  strength  parameters  of  the  upper  layer. 

The  locations  where  the  lower  layer  causes  a jump  in  the  normal 
stresses  are  determined,  and  the  approximate  value  of  the  jump  is 
computed.  In  the  second  step  the  interface  stresses  are  computed 
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with  the  strength  parameters  of  the  lower  layer  assuming  the  same 
centerline  geometry  as  in  the  first  step.  The  interface  stress 
distribution  for  the  two-layer  soil  is  then  determined  by  shifting 
the  normal  stresses  computed  with  the  strength  parameters  of  the 
lower  layer  by  the  appropriate  amount  within  the  locations  of  the 
stress  jumps  determined  earlier.  The  load,  drawbar  pull,  and  tor- 
que are  determined  for  this  interface  stress  distribution.  The 
procedure  is  repeated,  if  necessary,  for  a new  centerline  geometry 
until  the  computed  load  is  within  the  tolerance  limits  for  the  in- 
put load.  The  computation  scheme  is  shown  diagrammatically  in 
Fig.  23.  A more  detailed  flow  diagram  is  given  in  the  appendix. 

An  optional  provision  in  the  scheme  allows  the  consideration 
of  the  equalization  effect  of  opposite  directed  shear  stresses 
(discussed  in  Section  IV) . A linearly  varying  value  of  the  inter- 
face friction  angle  is  superimposed  on  the  uniform  8 normally 
assumed  in  tire-soil  interaction. 

A typical  interface  normal  stress  distribution  computed  by 
this  method  for  two-layer  soils  is  shown  in  Fig.  24.  The  following 
limitations  apply  to  the  use  of  the  model: 

• The  model  applies  only  to  cases  where  the 
deformed  centerline  of  the  tire  is  com- 
pletely in  the  upper  layer 

• The  strength  conditions  are  such  that  both 
the  cohesion  and  the  friction  angle  of  the 
upper  layer  is  either  higher  or  lower  than 
that  of  the  lower  layer  (Case  1 and  Case  2) 

The  pull  performances  estimated  by  this  tire-layered  soil  model 
were  compared  with  experiments  performed  at  WES  (Ref.  13).  Because 
of  the  difficulties  involved  with  the  preparation  of  layered  soil 
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Fig.  23  Scheme  of  Computations  for  the  Tire-Soil  Model 
for  Two-Layer  Soils 
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beds  where  the  strength  of  each  layer  is  controlled,  the  experi- 
ments at  WES  were  performed  in  a small  soil  bin  with  2.50-4  tires. 
Even  though  great  care  was  taken  to  prepare  each  layer  with  the 
specified  strength,  small  deviations  in  the  strength  properties  in- 
evitably occurred,  and  all  test  data  had  to  be  normalized  to  allow 
a meaningful  evaluation  of  the  effects  of  layering.  For  that 
reason,  the  normalized  test  data,  instead  of  individual  test  data, 
were  used  in  the  comparison  with  performance  predictions  by  the 
model. 


In  the  experiments,  the  differences  in  pull  performances  due 
to  the  lower  layer  were  found  to  be  relatively  small.  To  predict 
such  small  differences,  it  was  necessary  to  duplicate  the  pull  per- 
formance of  so-called  "base  line"  tests  closely.  In  the  develop- 
ment of  the  tire-soil  model  (Ref.  7),  various  relationships  that 
define  the  model  parameters  (limit  pressure,  deflection  coefficients, 
slip  coefficients)  were  developed  from  the  analysis  of  tire  per- 
formance tests  conducted  on  tires  4.00-7  and  larger.  For  the 
small  2.50-4  tire  used  in  the  experiments  in  two-layer  soils,  the 
model  parameters  estimated  on  the  basis  of  these  relationships 
yielded  predictions  generally  acceptable,  but  not  accurate  enough 
for  use  in  analyzing  small  differences  in  pull  performances  ob- 
served in  the  small  scale  mobility  tests  on  layered  soils.  Thus 
it  was  desirable  to  make  adjustments  in  these  parameters  so  that 
the  base  line  performances  could  be  more  closely  duplicated.  An 
uncertainty  arose,  however,  with  the  shear-slip  coefficients  in 
the  equation 


-(j+j  )/K 

j o' ' \ 

tan  6 = tan  6 l 1 - e ) 

max\  / 


T = T 


max 


(16) 
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Since  various  combinations  of  the  constants  i and  K can 

J o 

yield  the  same  6 value,  it  was  not  possible  to  establish  firm 
values  for  these  constants  from  the  small  scale  mobility  tests  that 
were  performed  only  for  one  value  of  the  slip  j,  at  approximately 
j = 20  percent.  To  avoid  the  problems  associated  with  establishing 
the  values  of  j and  K for  the  small  scale  mobility  tests,  it 
was  found  sufficient  for  the  comparative  analysis  of  these  tests  to 
establish  the  degree  of  shear  mobilization  defined  as: 

t tan  6 

u = = (17) 

t tan  6 

max  max 


The  degree  of  shear  mobilization,  as  defined  by  Eq.  (17),  was 
established  for  the  base  line  tests  and  retained  in  the  performance 
simulation  of  two-layer  conditions. 


In  the  tests,  the  strength  properties  of  the  upper  and  lower 
layer  were  determined  by  cone  penetration  tests;  the  values  of  cone 
indices  are  shown  in  the  respective  columns.  In  the  tire-layered 
soil  model  the  strength  properties  of  each  layer  are  characterized 
by  its  Coulomb  strength  parameters.  For  the  purpose  of  comparing 
the  measured  and  predicted  tire  performances  the  cone  index  values 
have  been  converted  to  Coulomb  strength  parameters  by  the  relation- 
ships established  for  the  Buckshot  clay  in  Ref.  7.  These  relation- 
ships are  as  follows: 


c(lb/sq  ft)  = 11.5  Cl 
cp(deg)  = 0.25  Cl 


(18) 


The  results  of  the  small  scale  mobility  tests  performed  in 
Buckshot  clay  on  various  combinations  of  upper  and  lower  layer 
strength  are  summarized  in  Tables  1 and  2 together  with  the  results 
of  simulation  by  the  tire-layered  soil  model. 
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TABLE  1 TIRE  PERFORMANCE  IN  TWO-LAYER  SOILS  — SOFT  OVER  HARD  LAYER 


i — i 

r-H  CM 
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67.4  t V I 0.0625  0.24  V V ¥ 0.23 


TABLE  2 TIRE  PERFORMANCE  IN  TWO-LAYER  SOILS  — HARD  OVER  HARD  LAYER 


24.7 


Table  1 shows  the  test  results  for  Case  1,  soft  over  hard 
layer.  The  base  line  test  results  are  shown  as  for  infinite  top 
layer  thickness.  The  gain  in  pull  performances  due  to  a stronger 
lower  layer  was  minimal  in  all  these  tests.  The  simulation  by  the 
tire-layered  soil  showed  no  increase  in  the  pull  performance  in 
Group  A (tests  1 through  4)  and  Group  B (tests  5 through  8).  The 
reason  for  the  stronger  lower  layer  having  no  influence  on  the  pull 
coefficient  in  these  simulations  is  that  neither  the  forward  nor 
the  rear  slip  line  field  computed  in  the  model  intersected  the 
layer  interface  in  these  cases.  The  question  arises  whether  the 
actual  failure  zones,  contrary  to  the  simulation,  intersected  the 
layer  interface,  or  whether  the  increase  in  pull  coefficient  was 
due  to  some  interaction  in  that  portion  of  the  upper  layer  where 
the  limit  pressure  governs  the  stresses  and  the  soil  is  not  in 
failure  conditions.  An  explanation  for  the  tire-layered  soil  model 
not  showing  the  increase  in  pull  coefficient  obtained  in  the  tests 
is  that  the  actual  friction  angle  might  have  been  higher  than  that 
estimated  by  Eqs.  (18).  While  the  over-all  performance  is  not 
sensitive  to  small  deviations  in  the  friction  angle,  the  depth  of 
the  slip  line  field  changes  appreciably  with  the  friction  angle  and 
could  thereby  affect  the  simulation  of  the  interaction  in  two-layer 
soils.  These  questions  can  only  be  clarified  by  an  experimental 
program  conducted  for  this  purpose.  At  any  rate,  a modification  of 
the  model  on  account  of  the  marginal  increases  in  the  pull  coeffi- 
cients in  these  tests  does  not  appear  justified.  The  simulation  of 
the  pull  performance  in  Group  C (tests  9 through  11)  and  Group  D 
(tests  12  through  16)  is  satisfactory,  even  though  the  simulation 
was  not  sensitive  in  the  case  of  the  largest  top  layer  thickness  — 
probably  for  the  same  reasons  as  cited  for  Groups  A and  B. 
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Table  2 shows  the  test  results  for  Case  2,  hard  over  soft 
layer.  In  this  case  the  tire-layered  soil  model  tended  to  over- 
predict the  decrease  of  pull  coefficient  due  to  the  lower  strength 
of  the  bottom  layer.  A look  at  the  test  results,  for  example  in 
Group  D,  casts  doubts  whether  the  tests  or  the  simulation  gives 
the  right  answer.  Attention  may  be  called  to  the  significance  of 
shear  stresses  at  the  layer  interface.  The  studies  performed  in 
connection  with  the  development  of  methods  for  constructing  com- 
posite slip  line  fields  showed  that  shear  stresses  at  the  layer 
interface  strongly  influence  the  bearing  stresses,  especially  in 
Case  2.  Such  shear  stresses  are  generated  by  the  stress  field  in 
the  upper  layer,  but  they  may  also  be  generated  and  "locked  in"  in 
the  process  of  the  preparation  of  layered  beds.  Locked  in  shear 
stresses  may  be  an  explanation  for  the  insensitiveness  of  the  tests 
to  the  weak  bottom  layer.  An  experimental  program  with  normal  and 
shear  stress  measurements  at  both  the  tire-soil  and  layer  inter- 
faces could  clarify  the  problems. 


VI.  THE  TIRE -LAYERED  SOIL  MODEL  AND 
THE  PROBLEM  OF  SLIPPERINESS 


In  connection  with  the  development  of  methods  for  the  con- 
struction of  composite  slip  line  fields  (discussed  in  Section  IV), 
it  was  pointed  out  that  the  construction  of  such  slip  line  fields 
may  not  be  feasible  if  the  strength  of  the  upper  layer  is  low  com- 
pared to  that  of  the  lower  layer.  Although  slippery  conditions 
have  not  been  defined  in  terms  of  strength,  it  is  reasonable  to 
assume  that  they  can  be  identified  with  the  condition  when  the 
construction  of  composite  slip  line  fields,  as  proposed  in  Sec- 
tion IV,  is  not  feasible.  At  this  point  it  is  useful  to  recall 
that  these  slip  line  fields  represent  the  solution  of  the  differ- 
ential equations  of  plasticity  for  soils.  Along  the  slip  lines, 
the  Mohr-Coulomb  yield  criterion  as  well  as  the  differential  equa- 
tions of  equilibrium  are  satisfied.  In  plasticity  theory  such  solu- 
tions are  labeled  as  "lower  bound"  solutions  that  become  true  solu- 
tions if  the  slip  lines  are  compatible  with  certain  kinematic  condi- 
tions . 

When  stress  solutions  do  not  exist  for  the  given  boundary  con- 
ditions, as  under  slippery  conditions,  the  upper  bound  theorem  of 
the  theory  of  plasticity  may  be  used  for  the  estimation  of  failure 
stresses.  According  to  this  theorem,  slip  line  fields  are  con- 
structed on  the  basis  of  kinematics,  and  the  stresses  are  integrated 
along  the  given  slip  lines.  Such  procedure  yields  a so-called  upper 
bound  solution.  Figure  16  shows  such  a kinematic  slip  line  field 
for  the  simple  bearing  capacity  case.  If  the  friction  angle  of  the 
slippery  layer  is  assumed  as  zero,  then  the  bearing  stress  can  be 
calculated  from  the  following  formula  (Ref.  11): 
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(19) 


p = c{\  ir  + ~)  + m2c 

where 

B = width  of  footing 
h = thickness  of  upper  layer 
ra  = factor  0 < m < 1 

An  inspection  of  Eq.  (19)  indicates  that  with  decreasing 
thickness  of  the  upper  layer  the  interface  pressure  p tends  to 
go  to  infinity.  In  the  case  of  slipperiness,  the  slushy  upper 
layer  may  be  squeezed  out  from  beneath  the  tread  to  some  extent, 
but  never  completely,  since  that  would  require  an  infinitely  large 
pressure.  According  to  the  concept  of  tire-soil  interaction  there 
is  a limit  to  the  pressures  that  soil  can  develop  under  the  tire 
load.  One  may  assume  that  under  slippery  conditions  the  upper 
slushy  layer  would  be  squeezed  out  from  underneath  the  tire  tread 
as  long  as  the  pressure  is  smaller  than  that  indicated  by  Eq.  (19). 
The  pressure  will  rise  due  to  the  gradual  decreasing  of  the  thick- 
ness of  the  upper  layer  as  the  slushy  material  is  squeezed  out,  but 
will  stabilize  at  the  limit  pressure.  Thus,  under  slippery  condi- 
tions the  interface  normal  stress  may  be  assumed  as  p^,  and  the 
interface  itself  flat,  paralleling  the  lower  hard  layer.  The  pull 
coefficient  may  be  approximately  computed  as: 


where  s is  the  strength  of  the  slush.  The  pull  coefficient  pre- 
dicted by  the  tire-layered  soil  model  for  slippery  conditions  will 
be  approximately  that  computed  from  Eq.  (20). 

Equation  (20)  denotes  the  pull  coefficient  for  slippery  con- 
ditions when  there  is  no  traction  component  from  the  lower  hard 
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layer.  The  lower  hard  layer  can  be  engaged  in  traction  only  if 
is  high  enough  to  cause  failure  and  penetration  in  the  lower 
hard  layer.  This  explains  why  tires  with  high  inflation  pressures 
work  better  under  slippery  conditions  than  low  inflation  pressure 
tires.  Slippery  conditions  may  occur  for  a wide  range  of  bottom 
layer  strength;  therefore,  the  use  of  higher  inflation  pressure 
would  be  advantageous  only  for  a specific  range  of  bottom  layer 
strength. 
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VII.  FIELD  DETERMINATION  OF  PROPERTIES 
OF  LAYERED  SOILS 

To  account  for  layered  soil  conditions  in  the  evaluation  of 
mobility  over  a terrain  that  exhibits  such  conditions,  it  is  neces- 
sary to  determine  in  the  field  the  variation  of  the  strength  proper- 
ties of  soil  with  depth.  The  suitability  of  the  present  field  tech- 
niques for  this  purpose  has  been  reviewed  and  the  possibilities  of 
improvements  investigated.  The  results  of  this  review  are  summarized 
in  this  section. 

PLATE  SINKAGE  TEST 

Although  the  plate  sinkage  tests  in  their  present  form  are  used 
primarily  for  the  establishment  of  pres sure -sinkage  curves,  a sharp 
break  in  these  curves  indicates  the  point  at  which  plastic  yielding 
starts  in  the  soil.  The  total  load  on  the  plate  at  this  point  may 
be  computed  from  t he  failure  stresses  chat  produce  plastic  zones  in 
the  soil  . The  plate  sinkage  test  differs  from  the  simple  bearing 
capacity  case  analyzed  in  Section  IV  in  that  it  is  an  axially  sym- 
metric case.  The  effects  of  layering  in  the  two  dimensional  and 
axially  symmetric  case  are  essentially  the  same.  In  a plate  sinkage 
test  in  a two-layer  system  the  bottom  layer  affects  the  bearing 
stresses  as  soon  as  the  plastic  zone  intersects  the  lower  layer. 

Thus  there  is  no  distinct  change  in  the  pres sure -sinkage  curve  that 
would  be  indicative  of  the  strength  properties  of  the  lower  layer, 
but  rather  a continuous  transition  in  the  pressure  that  would  be 
exceedingly  difficult  to  evaluate.  Thus,  in  principle,  the  plate 
sinkage  test  is  unsuitable  for  the  evaluation  of  the  strength 
properties  of  discrete  individual  layers. 
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CONE  PENETRATION  TESTS 


Cone  penetration  tests  are  intended  to  measure  strength  proper- 
ties of  soil  for  mobility  evaluation.  In  layered  soils,  not  only 
the  plastic  zones  created  by  the  cone  but  the  cone  itself  intersects 
the  layer  interface,  thereby  making  the  distinction  of  strength 
properties  of  individual  layers  impractical.  A few  cone  penetration 
tests  were  conducted  in  the  laboratory  in  layered  soils  to  demon- 
strate the  problem  and  see  whether  the  use  of  a right  angle,  in- 
stead of  a 30°  cone,  would  alleviate  the  situation.  Typical  re- 
sults of  these  tests  are  shown  in  Figs,  25,  26,  and  27.  Figure  25 
shows  the  variation  of  cone  penetration  with  depths  for  a 30°  and 
a 90°  cone  in  two  layers  of  Buckshot  clay.  With  the  30°  cone 
the  reading  for  the  1-in.  (2.5  cm)  thick  top  layer  (Cl  = 11) 

is  completely  obscured  by  the  presence  of  the  stronger  bottom  layer. 
This  is  understandable  since  the  more  than  3 cm  long  cone  pene- 
trates the  bottom  layer  even  when  its  base  is  at  the  surface.  The 
90°  cone  showed  momentarily  a reading  corresponding  to  the  strength 
of  the  upper  layer.  Nevertheless,  from  the  signature  itself  one 
could  not  evaluate  the  strength  of  the  upper  layer  even  if  the  thick- 
ness of  the  layer  were  known.  Figure  26  shows  the  variation  of  the 
cone  penetration  resistance  in  a Buckshot  clay  with  alternating 
layers  of  different  strengths.  The  changes  in  cone  penetration  re- 
sistance indicate  that  layering  is  present,  but  the  signature  in 
itself  is  insufficient  to  establish  depths  of  layer  interfaces  or 
the  strength  properties  of  individual  layers. 

Figure  27  shows  the  results  of  cone  penetration  tests  in  a two- 
layer  sand.  The  change  in  the  gradient  of  the  resistance  appears 
to  occur  at  the  layer  interface.  Another  change  in  the  gradient, 
at  point  "A",  may  be  the  result  of  the  slip  zones  reaching  the 
upper  loose  layer,  as  shown  by  dashed  lines  in  the  figure. 
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Pig.  26  Results  of  Cone  Penetration  Tests  In  Multilayered 
Clay  Soils 
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Fig.  27  Results  of  Cone  Penetration  Tests  in  Two-Layer  Sand 
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RING  SHEAR  TESTS 


From  the  preceding  discussions  it  is  obvious  that  field  tech- 
niques that  engage  both  layers  either  directly  or  indirectly  through 
the  creation  of  plastic  zones  at  some  distance  from  the  layer  inter- 
face are  unsuitable  for  the  measurement  of  strength  properties  of 
individual  layers.  The  ring  shear  test  that  has  been  used  in  mo- 
bility technology  for  some  tinE  engages  only  a horizon  of  the  soil 
fully;  it  therefore  appears  better  suited  for  determining  strengths 
of  individual  layers  than  do  the  other  techniques. 

The  drawback  of  the  ring  shear  test  technique  is  that  the  peak 
shear  stress  measured  at  an  applied  normal  stress  equals  the  shear 
strength  of  the  material  only  if  the  base  of  the  annulus  is  the 
failure  plane.  Observations  and  theoretical  considerations  indicate 
that  in  many  instances  this  condition  is  not  met.  The  failure  plane 
in  soils  is  perpendicular  to  the  plane  of  the  major  and  minor  prin- 
cipal stress.  If  the  lateral  stresses  in  the  soil  are  low,  then  the 
minor  principal  stress  is  in  the  lateral  direction  and  the  major 
principal  stress  is  in  the  plane  tangent  to  the  rotation  and  in- 
clined to  the  vertical  due  to  the  application  of  shear  stresses. 

The  failure  plane  is  an  oblique  one,  perpendicular  to  the  plane  that 
includes  the  minor  principal  stress  inclined  to  the  vertical.  The 
development  of  such  an  oblique  failure  plane  manifests  itself  by  an 
excessive  sinkage  or  "digging  in"  of  the  ring  shear  apparatus, 
since  plastic  flow  from  along  the  inclined  failure  plane  moves  soil 
out  from  underneath  the  annulus.  A detailed  discussion  of  this 
problem  is  given  by  Liston  in  Refs.  14  and  15.  He  proposes  an 
evaluation  procedure  that  takes  failure  along  the  oblique  plane 
into  account. 

The  development  of  an  oblique  failure  plane  nullifies  the  ad- 
vantage of  the  ring  shear  test  in  that  it  engages  only  a horizon 
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of  the  soil.  Thus,  for  the  purpose  of  measuring  the  shear  strengths 
of  individual  layers,  failure  along  oblique  planes  would  have  to  be 
prevented  rather  than  the  evaluation  procedures  modified  for  this 
condition.  Theoretical  considerations  indicate  that  failure  along 
oblique  planes  could  be  prevented  by  placing  a surcharge  (metal 
plates)  around  the  shear  ring  and  possibly  in  the  open  area  in  the 
center.  Preliminary  calculations  indicate  that  a surcharge  equal 
to  the  vertical  stresses  applied  to  the  shear  ring  would  be  suffi- 
cient for  this  purpose.  The  schematic  arrangement  of  a ring  shear 
test  with  surcharge  is  shown  in  Fig.  28. 

The  application  of  surcharge  would  prevent  any  appreciable 
sinking  of  the  ring  shear  device  during  a test.  To  reach  deeper 
horizons,  a cylindrical  bottom  split  scooper  would  have  to  be  at- 
tached to  the  device  to  remove  the  soil  to  the  desired  depth.  The 
ring  shear  test  would  then  be  carried  out  at  this  depth  to  deter- 
mine the  strength  properties  of  soil  at  this  horizon. 

Ring  shear  devices,  portable  and  track  mounted,  are  already 
available.  The  concept  of  surcharge  application  could  be  imple- 
mented with  minor  modifications  to  the  existing  equipment.  It  is 
strongly  recommended  that  a program  be  initiated  for  detailed  de- 
velopment of  the  modified  techniques  of  the  application  of  ring 
shear  devices  and  the  performance  of  laboratory  tests  for  the 
validation  of  the  concept.  The  development  of  this  improved  tech- 
nique would  solve  not  only  the  problem  of  determining  the  strength 
properties  of  layered  soils  but  could  also  be  used  to  determine 
strength  properties  of  cohesive-frictioned  soils  and,  therefore, 
fill  a gap  in  the  area  of  applicability  of  cone  penetration  tests. 
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NORMAL  STRESS  DUE 


Fig.  28  Schematic  Arrangement  of  Ring  Shear  Test  with  Surcharge 


VIII.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  variation  of  strength  properties  of  soil  with  depth,  fre- 
quently occurring  in  the  field,  and  its  effect  on  tire-soil  inter- 
action have  been  investigated  theoretically.  One  method  for  the 
consideration  of  a continuous  variation  of  strength  and  another  for 
the  consideration  of  discrete  layering  have  been  developed  and  ap- 
plied in  an  expanded  tire-soil  model. 

In  the  tire-soil  model  expanded  for  continuously  varying  soil 
strength  a linear  variation  of  soil  strength  with  depth  is  con- 
sidered. Analyses  performed  by  this  model  indicate  that  the  pull 
performance  depends  more  on  the  strength  of  the  surface  than  at 
depth  as  long  as  the  rate  of  decrease  of  soil  strength  is  less 
than  a certain  limit.  At  this  limit  a no  go  situation  occurs  quite 
abruptly.  It  is  recommended  that  the  "constant  strength"  concept 
requiring  an  averaging  of  the  soil  properties  within  a certain 
depth  be  changed  to  a "linearly  varying  strength"  concept  and  that 
a constant  plus  a rate  of  variation  be  used  for  the  characteriza- 
tion of  soils  in  conjunction  with  the  variable  strength  tire-soil 
model.  The  accuracy  of  performance  prediction  could  be  signifi- 
cantly improved  thereby  with  little  additional  effort  in  field 
data  processing. 

The  consideration  of  discrete  layer  interfaces  in  plasticity 
theory  solutions  required  both  extensive  theoretical  research  and 
the  development  of  new  numerical  techniques.  These  applied  to  the 
simple  bearing  capacity  case  resulted  in  a new  method  of  construct- 
ing composite  slip  line  fields  and  obtaining  bearing  stresses  for 
two-layer  soils.  An  approximate  method  for  estimating  these  stresses 
was  applied  in  the  tire-soil  model  for  the  consideration  of  two- 
layer  soils.  The  simulation  by  the  model  of  small  scale  tire  tests 


performed  in  two-layer  soils  was  acceptable,  but  some  assumptions 
in  this  two-layer  model  require  further  experimental  research.  It 
is  recommended  that  an  experimental  program  be  initiated  in  which 
both  tire-soil  and  layer  interface  stresses  are  measured „ Such  a 
program  would  be  not  only  useful  for  the  validation  of  the  concepts 
adopted  in  the  tire-layered  soils  model  but  would  also  contribute 
to  the  understanding  of  the  basic  phenomena  in  tire-soil  interac- 
tion in  layered  soils. 

The  techniques  employed  for  the  determination  of  soil  proper- 
ties in  the  field  were  examined  with  respect  to  their  suitability 
for  distinguishing  the  strength  of  individual  layers.  Neither  the 
plate  sinkage  nor  the  cone  penetration  tests  are  suitable  for  this 
purpose.  The  ring  shear  test  is  best  suited,  but  the  application 
of  a surcharge  is  required  to  prevent  lateral  failure.  It  is 
strongly  recommended  that  the  existing  ring  shear  devices  be  adapted 
for  the  accommodation  of  surcharge  plates  and  the  procedure  be 
validated  in  an  experimental  program.  The  proposed  procedure  is 
expected  to  be  applicable  not  only  in  layered  soil  conditions  but 
also  for  determining  strength  properties  of  frictional-cohesive 
soils  for  which  present  procedures  are  inapplicable. 
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APPENDIX 

FLOWCHARTS 
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FLOWCHART  FOR  THE  COMPUTATION  OF  SLIP  LINE  FIELDS 
FOR  SOIL  STRENGTH  VARYING  WITH  DEPTH 


69 


10/30/75 


AUTOFLOW  CHRRT  SET  - FLOCHRRT  SLF00010 


PAGE  01 


CHRRT  TITLE  - SUBROUTINE  SLFt 


Z_J£ H7 


SUBROUTINE  SLIP 
COMPUTES  COORDINATES, 
SIGMA  AND  THETA 
VALUES  OF 

SLIP  LINE  FIELDS  FOR 
INPUT  DATA 

TRH NS HITTER  FROM  THE 
MAIN  PROGRAM. 

IT  RETURNS  TO  THE 
HA IN  FRO GRAM  CENTRAL 
ANGLES  (HUH, 

X COORDINATES  l U IJM  , 
Z COORD  I ANTES  IV  [JJJ 
OF  POINTS  AT  THE 
INTERFACE  RNu  THE 
COMPUTED 
NORMAL  STRESSES 
(OM'i  ) AND  SHEAR 
STRESSES  (E  I'J)  RT 
THESE  POINTS 
INITIALIZE  ARFiATS 


71 


1 D/30/75 


nUTOFLOW  CHART  SET  - FLOCHRRT  SLFOOOIO 


PAGE  02 


CHART  TITLE  - SUBROUTINE  SLFI 


“02. 15* 


— >■ 


/5LF012goXXv18 
/ 

/c 


FI  .LE.  0.0 


1 7 1 


F1MAX  = 

AH  AX  l (F 1MRX,  FI) 


F1MIN  = 

AM  IN  l (FIMIN,  FI) 


.TRUE 


Tfaise 


/3LF0123CI/ 


19 


1 a?. 

HR.'TE  TO  DEV  / 

3 / 
VI H FORHRT  / 
SOOT  / 
FROM  THE  LIST  / 


/SLF03.310/ 


NOTE  2T 


LIST  = FI,  Cl. 
XX.  I.  J.  ZiI.J) 

■~T' 


3£0/_[ 25_ 


/SIF05320/  . 

/ WRITE  TO  DEV  ,7 
/ 10  / 

/ VIA  FORMAT  / 

/ 242  / 

/ FROM  THE  LIST  / 


/SLf0;320,' 


MOTE  26 


LIST  = XX,  I,  J, 
H a,  Ji 


7 


27 


02.  c'0>“ — ’> 

™ J ^ 

/ HR!  TL  TO  DEV  / 
/ 3 / 

/ VIA  TOCNAT  / 

/ Si;  55  / 

/ FR!1M_TH£  LIST  J 


/SLFOTSSO/  ! 

Lm£_23 

~ - — -I 

[|  xy.  [>  j. 

/SLF03370' 

2 _9_ 

/ WRITE  TO  DEV  / 
/ 1(1  / 
/ VIA  FOF  HOT  / 

/ / 

/ FF;(1M  THF  1 1ST  / 


Z 

/SLF0537 


NOTE  SO 


[LIST  = XX,  I.  J,  n 
Z(IpJ)  

250 

31 

I H I = 0 

0*i . I 3*r  — i* 
1 00  ci 

( EXIT  ) 


72 


10/50/75 


flUTOFLBW  CHRRT  SET  - FLOCHRRT  SLFOOOIO 


PfiGE  03 


CHART  TITLE  - SUBROUTINE  SLH 


COMPUTE  CO  8 RO  MATES 
WITH  FIRST  ESTHETES 


73 


10/30/75 


RUT0FLPW  CHRFIT  SET  - FLOCHRRT  SLFOOOIO 


PR9E  01 


CHfiRT  TITLE  - SUBROUTINE  SLF[ 


COMPUTE  COOROIHRTES 
AND  VARIABLES  AT 


PROGRAM 


10/30/75 


RUTOFLOH  CHRRT  SET  - FLOCHRRT  SLFOOOIO 


PAGE  05 


CHARI  TITLE  - SUBROUTINE  SLF1 


75 


10/30/75 


RLFTOFLDW  CHRRT  SET  - FLOCHRRT  SLFOOOIO 


PAGE  06 


CHART  TITLE  - SUBROUTINE  SLFt 


10/30/75 


RUTOFLnH  CHART  SET  - FL0CHRRT  SLFOOOIO 


PRGE 


07 


CHART  TITLE  - SUBROUTINE  SLF[ 


77 


JO/30/75 


AUTOPLDW  fHRRT  SET  - FLOCHflRT  SLFOOOIO 


CHART  TITLE  - SUBROUTINE  SLFI 


’07. 03k- 
170 


HI  = HfJ'i  - H C J 

n 


H2  = BH  - H (J  - 

I) 


H3  * H2/HI 


HU  = H3w  fX  Tf-J  + 

1 - J.JJ  - X IN  - 
J ♦ 2,J  - 1 1 ) 

XfH  * 1 - J,J)  > 
X IN  - J + 2..J  - 
1)  HU 


/SI.FCS306-0/"  \.  12 


FLOWCHART  FOR  COMPUTATION  OF 
TIRE  PERFORMANCE  IN  TWO-LAYER  SOILS 


79 


1 J/D5/75 


flUTOFLDH  CHART  SET  - FLOCHART  LESOOOIO 


PRGE  02 


CHART  TITLE  - PROCEDURES 


81 


1 3/D5/75 


flUTOFLOH  CHfiFiT  SET  - FLOCHRRT  LESOOOIO 


PAGE  03 


CHART  TITLE  - PROCEDURES 


02.21“ ?| 

/LE SO 0710/ 

/ READ  FROM  DEV 

/ 5 
/ VIA  FORMAT  j 
/ 1154D  / 

/ INTO  THE  LIST  / 


/ LE son 7 10/  I NOTE  02 

rnsTTsTliSr 


03.19 n 

555  03 


/ WRITE  TO  DEV 
6 

VIA  FORM ft T j 

SOI  / 


/WRITE JO  OEV 

/ VIA  FORMAT  / 
/ '13?0 


/LESO  0 760/'  05 

/ READ  FROM  DEV 

/ VIA  FORMAT  / 

/ 'V.M  / 

/ IMTO  THE  I.JST / 

/LESO 0780/  NOTE  06 


jTnf  = co:.  foT  rpj 


/ WRITE  TO  DEV 

/ 6 , 

/ VIA  FORMAT  / 

502  / 


/WRITE  TO  DEV 
/ -2 
' VIA  FORMAT  / 
1575  / 


/ READ  FROM  OEV 

/ VIA  FORMAT  / 
/ 4540  / 

£ INTO  THE  LIST  / 


/LESO OS 10/  1 NOTE  10 


LIST  * C02,  F0£\ 
002 


DETERMINE  RELATIVE 
STRENGTH  OF  LAYERS* 
UFTE  LATER  MEAKEFi 


82 


1 J/D5/75 


fiUTOFLBH  CHRRT  SET  - FLOCHRRT  LESOOOIO 


PAGE  04 


CHRRf  TITLE  - PROCEDURES 


83 


1 J/05/75 


flUTOFLBW  CHRRT  SET  - FLOCHRRT  LESOOOIO 


PAGE  05 


CHART  TITLE  - PROCEDURES 


1 J/D5/75 


flUTOFLBH  CHRRT  SET  - FL0CHRRT  LESOOOIO 


PRGE  06 


CHRRf  TITLE  - PROCEDURES 


85 


11/05/15 


AUTOFLDH  CHRRT  SET  - FL0CHRRT  LESOOOIO 


PAGE  01 


CHART  TITLE  - PROCEDURES 


/LE5021HD/ 

XXX  = 

- 1 

DOF  = 

OF 

ZDEP 

zLfir 

JPAS7 

= rr 

ZPftSI 

- 0 

ZLflSl 

* 0 

/LE50220D/ 


OH 


SLFI 


/LE50221D/^x  D5 


I HI  .EQ.  0 


'\TP.IJE 


FALSE  | 22 
1H 


/IE 502220/ 

QD  « PPL  - OOF 
Ml  = JJ1 


39  3 
06 


/LES022Ujl/-  \ 07 


[LAY  .EO.  2 
'l  FALSE 


\ TRUE 


60 


CHECK  FOR  HATCH  OF 
NORHRL  STRESS  AT 
SEPARATION  ANCLE 
irjQFl 

AND  LIHIT  PRE5SURE 
FPL  ) 


/LE502270/  X OR 
/ X'\ 

TRUE,. 


N. 


QO  .LT.  0 

[false 


8.01 
-^7 


86 


/D5/75 

HnRF  TITLE  - PROCEDURES 


nUTOFLOM  CHRRT  SET  - FLOC'HRRT  LESOOOIO 


PAGE  08 


11/D5/75 


RUTOFLOW  CHRRT  SET  - FLOCHRRT  LESOOOIO 


PRGE  10 


CHRRT  TITLE  - PROCEDURES 


IF  NORHRL  STRESS  FDR 
MAX.  ENTRT  ANGLE  LESS 
THAN  LIMIT  PRESSURE 
MAKE  LIMIT  PRESSURE 
EQUAL  TO  THAT  NORMAL 
STRESS  RAD  REFER! 
COMP. 


09.02 ? 

79 

L 01 

PPL  = QQF  | 

07.07 * 7> 

80 

02 

V9  = 

VV  IJ-Jl) 

U9  - 

UU  IJJ1) 

X9 

* VVK 

29 

» UUK 

N2  = 

37  - N1 

STORE  VALUES  OF 
INTERFACE  STRESSES 
RND  GEOMETRT 


BEGIN  DO  LOOP 
110  N = N2.  36 

on 

HHH  (N)  = HH  (37  - 
N) 

uGQ  (N)  = Qfl  (37  - 

m 

EEE  CN)  = EE  (37  - 
N) 

WV  (N)  = VV  (37  - 
N) 

/LE5027W 

OS 

ULIU(N)  = UU  (37  - 
N) 

no/1 

NO  ,/ 

./ 

/ 

END  OF  OD 
LOOP?  , 

\ / 
/ 
yes 


FALSE 

/LE50279O/ 

09 

AAOFI 

= RAO 

flflMFI 

= RAM 

COMPUTE  ANGLE  OF 
INTERFERENCE  fHLRTlJ 
AND  NORMAL  STRESS 
THERE- 

(QLA'fl)  FOR  FRONT 
FIELD 


/LES0282O/ 


LO 


2RHT  = IZF'AST 
ZDEP) / (ZPAST  - 
ZLAST) 

HLAT1  = 
HHCJPASTI  - 
(HH  UP  AST)  - 
HHfJPAST  - 

1)  :i«zrat 


/LESD28U0Z 


DLRY1  = 

QQ  (JPHSTI  - 
(QGCJPAST)  - 
QQUPAST  - 
1) ) *ZRAT 

SlGLfll  = SIGLA 


12 


JLRT  = 0 


10.07*- 

89  I NOTE  13 


BEGIN  DO  LOOP 
1)2  N = 1.  HI 


H 


hh  m 

= 0 

QQ  CN) 

= 0 

EE  CN) 

= 0 

VV  (N) 

= 0 

UU  (Ml 

= 0 

“V  END  OF  DO  ^ 
X LOOP? 


YES 


0 

^^120 


89 


11/D5/75 


RUTOFLOW  CHART  SET  - FLOCHART  LESOOOIO 


PRGE  11 


CHART  TITLE  - PROCEDURES 


15/05/75 


RUTOFLDW  CHART  SET  - FLOCHRRT  LESOOOIO 


PRGE  12 


CHRRT  TITLE  - PROCEDURES 


if 


/LES03580/ 

7 

HRAT  = (ZPRST  - 
ZDEP)  / (ZPRST  - 
ZLRST) 

HLRT2  = 
HHIJPA3T)  - 
(HH  C JPflST}  - 
HH  IJPAST  - 
m «ZRAT 

/LE SO 3 600/ 

02 

QLAY2  * 

QQ  IJPAST)  - 
(QQ  UPftST)  - 
QQ  IJPAST  - 

m hzrat 

SIGLR2  = S1GLR 


HLAY2  = - HLAY2 


11.18- 
/LE50346D/  I 


| XXX  = 3 | 

11,18  1 

m 

ZDEP  = ILRY  - 
(29  - ROxCSSIftun 

11.17 7 

117 

05 

DDF  = DF 
JPflST  = 0 
ZPABT  = 0 
ZLRST  = D 


/LES0352QXX  08 


TRUE 


/LES0355DXX  OR 

X \ 

X \TRUE 

\HPflST  . EQ.  0 > 

\X 


COMPUTE  ANGLE  OF 
INTERFERENCE  IN  REAR 
FIELD  (HLOT21 
AND  NORMAL  STRESS 
THERE  (QLAY5) 


1 


ii 


7 


WRITE  70  DEV  / 
6 

HR  FORMRT 
SIS 


0 


99 


131 


12 


KLflY  = 0 


12.02  — 

133  /X  13 


<^XXX  .GE.  2 

\x 

FALSE 


\ 14 


X 


QQF  .LT. 
V I.  95*PPU 


TRUE 


FRL5E  (13 


121 


15 


.TRUE 

QQF  ,LE.  PPL 

\X 

\X  . . 

FOLSE  f 13  ^ 
^09 J 

130 


LIMIT  NORMAL  STRESSES 
TO  i HE  VALUE  OF  LIMIT 
PRESSURE 


/LES03700/  HATE  IS 


FAL5E/  V 

K QQ  IN)  . GE.  > 

X FPL  / 

XX 

TRUE 


/LE SO 3 720/ 


19 


QQ  IN)  = PPL 


NOTE  20 


flBS  (EE  (N) ) .GT. 

( (QQ  (N)  + 

PS  1 1 xTflN(Fin 


/LES03730XX  21 

_X  \, 

X 


v 


.TRUE 

v SEE  NOTE  X 
X ABOVE  X 

\ X" 

XX  r 

IfALSE  ( 13  \ 


U/D5/75 


flUTOFLDM  CHRRT  SET  - ELOCHRRT  LESOOOIO 


PAGE  H 


CHRRT  TITLE  - PROCEDURES 


13.16- 
/LESOUllOy  X 03 


TRUE 


/LES0U12D/ 


| N31  = N3 

/LES0U13D/ 

NOTE  05 

BEGIN  DO  LOOP 
161  N = N3.  36 

08 

HA 

on  = 

HHH 

INI 

QA 

on  = 

DQQ 

(HI 

ER 

(N!  = 

EEE 

!H( 

VN 

on  = 

YVV 

INI 

/LESGUIsn/ 

07 

I UN 

IN)  = 

UUU 

INI 

SET  LOWER  LAYER 
STRENGTH  RND  COMPUTE 
INTERFACE 
STRESSES  THEREFOR 


/LESCHi220/''\  09 

/ \ 

true/  X 

n JLflY  .EQ.  0 > 


v 

I FALSE 


BEGIN  DO  LOOF 
163  N = N3.  36 

7 

163  / 
NO  / 

.X 

\ LOOP? 


/ 


YES 

/LES0U38O/ 

12 

P3IF  = 

C01/IAN  IFCUmHOI 

PSIR  = 

COR  1 /TAN  (FOR  Ik  HO) 

NS  2 = 

N3  | 

CONFUTE  FINAL 
INTERFACE  STRESSES 
FOR  TWO-LATER  SYSTEM 

/LE 30 4400/  NOTE  13 


BEGIN  DO  LOOF 
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IMPLICIT  REAL  (L) 

COMMON  RO, GO, CO, FI .01 ,R0, RM, fll , VK,UK,PL, DF. XX, 

QF.RF.VA, UR,  J1 , H (36)  ,Q  136)  ,E  (36)  ,11(36)  ,1  (36) 

COMMON  ZOEP , ZPA3T, ZLAST , JPAST, ISTREN.IHI 
COMMON  F02, F0R2, C02, CORO. HO, SIGLR.KDET 
EQUIVALENCE  (RRO.RO)  , (GOO,  GO)  , (CCO,  CO)  , (FF1  ,FU  , 

(001,01),  (RRO.AO)  , (RAM, AM),  Iflfll.Rl),  (VVK.VK)  , 

(UUK, UK) , (PPL,  PL)  , (DDF, OF) , IXXX.XX) . (OQF.QF)  , 

(RAF.RF),  (WR.Vfll  . (UUfl.UR),  IJJ1.J1),  (HH(1)  ,H(!)!, 

(ClO  (1) , 0 (I) ) , (EE  (1)  , E (1) ) , (UU  (H  ,U  (1)  I , (W  ID  , V II) ) 

DIMENSION  QQ!36)  ,EE(36)  ,HH  (36) 

DIMENSION  UU  (36)  , VV  (361 
DIMENSION  CIQQ  (36) , EEE  (36)  ,HHH(36) 

DIMENSION  UUU  (36' . VVV  (36) 

DIMENSION  OR  (36)  .Efi  (36)  .HR  (36) 

DIMENSION  VN  (361  .UN  (36) 

STATEMENT  FUNCTION  DEFINITION:  CELT  (D.F)  =RTHN  (SIN  (SIN  (0)  /SIN  (FI)  /SORT  (1-  (SIN(D)  /SIN  IF))  **2)) 

STATEMENT  FUNCTION  DEFINITION:  QUA  (09 ,F9)  = COS  (D9)  »SQRT  (COS  (D9)  **?-C0S (F91  x«2) 

STATEMENT  FUNCTION  DEFINITION:  CHIP  (D9.F3)  =035  (D9)  -SORT  (COS  (09)  mh2-C0S (F9)x«2) 

STATEMENT  FUNCTION  DEFINITION:  EPO  (F9.T9.T8)  -EXP (2x  (TS-T8)  *S IN  IF9) /COS  (F9)  I 

STATEMENT  FUNCTION  DEFINITION:  TAN  (F9)  =SIN  IF9)  /COS  (F9) 

STATEMENT  FUNCTION  DEFINITION:  TANO  (F9j. >SIN  (2*PI*  (F91  /36CI)  /COS  (2xPIx  (F9) /360) 

STATEMENT  FUNCTION  DEFINITION:  RET  (ES.D9.A8,  ft 9)  -A10G  (E3)  / (09-A8-R9) 

3 FORMAT  (1H0.  ’TIME*  ’,[61 

503  FORMAT (1H  INPUT  RADIUS, WIDTH  IN  FT  ’ I 

4545  FORMAT  (’  RA.BCi’) 

45(40  FORMAT  (3F  10. 31 

SOU  FORMAT  (1H  INPUT  LOAD  (LBS).  P3UBL  (PSD  ' ) 

4550  FORMAT!'  LD.PO'I 

505  FORMAT  Uh  , 'INPUT  DEFL  .COEFF, DELTA  COEFF  ’ ) 

4555  FORMAT  (’  DE.DF'I 

506  FORMAT  (1H  ,'  INPUT  SLIP  PARAMETERS  JZEP.O  AND  K ’) 

4580  FORMAT!'  SJ.SK’I 

SO?  FORMAT (1H  , ’ INPUT  SLIP.  DEPTH  OF  SECOND  LATER’  ) 

4565  FORMAT!’  SL.ZLAV’) 

501  FORMATOH  . 'INPUT  UPPER  LATER  COHESION.  FR.  ANGLE. GAMMA  ') 

4570  FORMAT C COl ,FOI , GO  1 ’ ) 

502  FORMAT (1H  .'  INPUT  LONER  LATER  COHESI ON , FR . ANGLE, GAMMA  ') 

4575  FORMAT  ( ’ C02 . FfJ2 , G02  ' ) 

511  FORMAT ( 1H  .’  RELATIVE  STRENGTH  OF  LATER  INCONSISTENT, TRY  AGAIN’) 

12  FORMAT (1H  , 'RAD. = ’. F9. 2. 1 WIDTH*’ .FS.2. ’ LOfiD= ’ , F9. 2. ’ 

PSUBL-  ’.F9.2) 

IN  FORMAT (1H  , ’DEFL.  COEF.= ’,  F6.  3.  ’ DELTA  COEFF.  =’,F8.3  ) 

16  FORMAT (1H  , ’DEPTH  OF  2M0  LATER-  ’, FI  0.31 

19  FORMAT ( 1H  . ’J  ZERO  = ’.F8.4,’  K = ’.F8.U.’  SLIP  = ’.FS.U) 

18  FORMAT  (1H  ,’  DELTA  = ’.F9.3,’  DEG.  OF  SHEAR  MOB=  ’,F9.3) 

13  FORMAT ( 1H  , ’ 1ST  L.  COH=  ’ , F 1 0 , 3 , ’ PH  I = ’,F8.3, ’ GAMMA-  ’,F8.2) 

16  FORMAT  (1H  , ’2ND  L.  COH=  ’,  F 10.3,  ’ PH  I = \FtV3,’  GAMMfl=  ’ , FS . 2) 
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SUBROUTINE  SLFI 
COMPUTES  COORDINATES, 
SIGMA  AND  THETA 
VALUES 

FOR  A SLIP  LINE  FIELD 
FOR  GIVEN 
ANGLES, INTERFACE 
FRICTION 
ANGLE  AND  LIMIT 
FRESSURE.  THE 
COMPUTED  INT.  FREE 
NORMAL  STRESSES 
(D(J)I,  SHEAR 
STRESSES 

(EIJ)  ) .COORDINATES  OF 
(NT.  FACE  (VUi.UUO 
PERTAINING  TO  ANGLES 
HtJ)  ARE  TRANSMITTED 
TO  THE  HAIN  PROGRAM 
THE  INTERSECTION  OF 
THE  SLIP  LINE  FIELD 
WITH  THE  INTEFILATER 
BOUNDARY  AND  THE 
SIGMA  VALUE  FOR  THE 
HAT.  INT.  FACE  SHEAR 
STRESS 

(SIGLR)  IS  ALSO 
DETERMINED 
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DISTRIBUTION  LIST 
(As  of  March  1975) 


DRdt/I- 

Please  notify  USATACOM,  -AM5TA-RHM,  Warren,  Michigan  48090, 
of  corrections  and/or  changes  in  address. 

No.  of 
Copies 

50  Commander 

U.S.  Army  Tank-Automotive  Command 
Warren,  MI  48090 

ATTENTION: 

Chief  Scientist,  AMSTA-CL  (230) 

Director  of  RD&E,  AMSTA-R  (200) 

Deputy  Director  of  RD&E,  AMSTA-R  (200) 

Foreign  Intelligence  Ofc,  AMSTA-RI  (200) 
Systems  Development  Division,  AMSTA-RE  (200) 
System  Engineering  Division,  AMSTA-RB  (200) 
Engineering  Science  Division,  AMSTA-RH  (200) 
Armor  & Components  Division,  AMSTA-RK  (215) 
Methodology  Function,  AMSTA-RHM  (215) 
Propulsion  Systems  Division,  AMSTA-RG  (212) 
Canadian  Forces  Liaison  Office,  CDLS-D  (200) 

US  Marine  Corps  Liaison  Ofc,  USMC-LN0  (200) 
Project  Manager  MICV,  AMCPM-MCV  (Dqr) 

Product  Manager  XM861,  AMCPM-CT  (Dqr) 

Project  Manager,  M60  Tank  Dev  (AMCPP-M60TD)Dqr 
Project  Manager  XM-1 , AMCPM-GCM  (dqr) 

Technical  Library,  AMSTA-RWL  (200) 

TRAD0C  Liaison  Ofc,  TRADOC-LNO  (200) 

2 Chief,  Research  Development  and  Acquisition 

Department  of  the  Army 
Washibgton,  DC  20310 

1 Superintendent 

US  Military  Academy 
ATTN:  Professor  of  Ordnance 

West  Point,  NY  10996 
>• 

1 Commander 

US  Army  Logistic  Center 
ATTN:  ATCL-CC  (Mr.  J.  McClure) 

Ft.  Lee,  VA  23801 


1 


1 


2 


3 


3 


3 


2 


1 


Commander 

US  Army  Concept  Analysis  Agency 
Long  Range  Studies 
8120  Woodmont  Avenue 
Bethesda,  MD  20014 


Oept  of  the  Army 
Office  Chief  of  Engineers 
Chief  Military  Programs  Team 
Rsch  & Dev  Office 
ATTN:  DAEM-RDM 

Washington,  DC  20314 

Director 

US  Army  Corps  of  Engineers 
Waterways  Experiment  Station 
P.0.  Box  631 
Vicksburg,  MS  39180 

Di rector 

US  Army  Corps  of  Engineers 

Waterways  Experiment  Station 

ATTN:  Mobility  & Environmental  Laboratory 

P.0.  Box  631 

Vicksburg,  MS  39180 


Director 

US  Army  Cold  Regions  Research  & Engineering  Lab 
ATTN:  Dr.  Freitag,  Dr.  W.  Harrison,  Library 

P.0.  Box  282  ^ 

Hanover,  NH  03755  ^ 


Commander 

US  Army  Materiel  Command 
AMC  Buildipgf,  Room  8S56 
ATTN:  Mr/R.  Navarin,  AMCRD-TV 

5001  Eisenhower  Avenue 
Alexandria,  VA  22333 


President 

Army  Armor  and  Engineer  Board 
Ft.  Knox,  KY  40121 


Commander 

US  Army  Arctic  Test  Center 
APO  409 

Seattle,  Washington  98733 


114 


2 Commander 

US  Army  Test  & Evaluation  Command 
ATTN:  AMSTE-BB  and  AMSTE-TA 
Aberdeen  Proving  Ground,  MD  21005 

2 Commander 

Rock  Island  Arsenal 
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